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ABSTRACT 

A sample of 196 AGNs and ULIRGs observed by the Infrared Spectrograph (IRS) on Spitzer is analyzed 
to study the distribution of the strength of the 9.7^m silicate feature. Average spectra are derived for quasars, 
Seyfert 1 and Seyfert 2 AGNs, and ULIRGs. We find that quasars are characterized by silicate features in 
emission and Seyfert Is equally by emission or weak absorption. Seyfert 2s are dominated by weak silicate 
absorption, and ULIRGs are characterized by strong silicate absorption (mean apparent optical depth about 
L5). Luminosity distributions show that luminosities at rest frame 5.5/xm are similar for the most luminous 
quasars and ULIRGs and are almost 10'' times more luminous than the least luminous AGN in the sample. 
The distributions of spectral characteristics and luminosities are compared to those of optically faint infrared 
sources at z 2 being discovered by the IRS, which are also characterized by strong silicate absorption. It 
is found that local ULIRGs are a similar population, although they have lower luminosities and somewhat 
stronger absorption compared to the high redshift sources. 

Subject headings: galaxies: active — galaxies: quasars — galaxies: ISM — infrared: galaxies 
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1. INTRODUCTION 

Previous observa tions with the Iri frared Spectrograph 
(IRS)"^ on Spitzer (iHouck et al.l l2004h have discovered a 
high redshift population of optically faint infrared sources. 
This population has been derived from sources having ex- 
treme infrared to optical flux ratios (IR/opt > 40, IR/opt = 
i^f,y(24/im)/i^f !,(/?)) as selected from sur veys at 24/xm with the 
Multiband Imaging Photometer (MIPS. Rieke et al.l l20b4) on 
Spitzer. Redshifts were determined from mid-infrared spec- 
tral features for 43 of 58 sources observed in the Bootes field 
diLouck et al. 2005; Weedman et al. 2006c) and for 14 of 18 
radio sources in the Sp itzer First Look Survey (FLS) field 
d Weedman et al.ll2006bh . For the 57 sources which have red- 
shifts, the redshift (z) ranges from 0.7 to 2.7, with a median 
z = 2.2. 53 of these 57 sources have silicate absorption, with 
only 4 having redshifts derived from polycyclic aromatic hy- 
drocarbon (PAH) emission. This is in contrast to spectral fea- 
tures of sources selected on the basis of their near-infrared or 
submillimeter characteristics of the spectral energ y distribu- 
tions, most of which show PAH emission features dYan et al.l 
12005; Lutz et al. 2005; Weedman et al. 2006aT). 

The presence of strong silicate absorption and the absence 
of PAH features in samples chosen from extreme values 
of IR/opt have been interpreted to mean that these sources 
are heavily obscured AGN. Extreme IR sources with con- 
firmed j)gtical_AGN_properties were also found by ISO- 
CAM (Haa s et aLll2004l) . It is therefore essential to compare 
the characteristics of these high redshift, heavily absorbed 
sources to known examples of AGNs and Ultraluminous In- 
frared Galaxies (ULIRGs) to determine if there are local ana- 
logues to the obscured, high redshift sources. In this letter, 
we examine this issue by studying the distribution of silicate 
strengths of local AGNs and ULIRGs. 
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The distribution of the mid-IR properties of different 
types of known AGNs has been discussed in various studies 
(Heckman et al. 1994; Haas et al. 2003; Siebenmorgen et al] 
2004; .Clavel et alJ 120001: ISM et al.1 120061: IBuchanan et all 
.2006) . The mid-IR spectral features, especially the sili- 
cate features at 9.7 and 18/im are partic ularly important i n 
the context of the unification model (eg. 'A ntonuccilll993h . 
The unambiguous detection of silicate emissions in quasars 
with Spitzer is a strong support for the un ification model 
(ISiebenmorgen et al.l l2005t lllao et al.l l2005h . but there is a 
large variety of silicate features i n AGNs. Silicate ab sorp- 
tion is seen in some type I AGNs ('Weedman et al."2005) and 
emission in some type II AGNs (Sturm et al. 2005). It is im- 
portant, therefore, to investigate the distribution of silicate 
strengths within a large sample of different types of AGNs. 
The distribution will constrain the models of AGNs (eg. 
'Efstathiou2006'; 'Fritz et al."2006'; 'Dullemond & van BemmeJ 
2005; Nenkova et al. 2002; Pier & KroHk 1993, 1992) and in 
turn provide insights into the geometry of the dusty structure, 
the optical depth of the clouds, and the filling factor along the 
line of sight. 

2. OBSERVATIONS AND DATA REDUCTION 

We gathered a large sample of local AGNs and ULIRGs 
from archival and published literature that have been ob- 
served with the Spitzer IRS. The AGN sample com- 
bines the unpublished GTO archival data (progra m 14) 
and pu blished data from Hao et al. (2005), W eedm an et 
(2005), IBuchanan et aH t2006) and IShi etalJ (l2006l) . The 
ULI RG sample is drawn from the GTO p rogram (program 
105).( lArmus et al.ll2007l: ISpoon et al.ll2007l) . The total sam- 
ple includes objects having a large variety of classifications: 
radio-quiet QSOs, radio-loud QSOs, Seyfert Is and Seyfert 
2s, IRAS-discover ed sources, and 2MASS-selected red AGNs 
(ISmith et alJl2002i) . 

We use only the IRS observations of sources which in- 
clude all low-resolution modules: Short-Low2 (SL2, 5.2- 
7.7/^m), Short-Lowl (SLl,7.4-14.5Aim), Long-Low2 (LL2, 
14-21. 3/im) and Long-Lowl (LLl, 19.5-38/im). We also ex- 
clude objects that have redshifts larger than 0.5, so that we 
can compare characteristics in similar rest frames. 
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To define sub samples of objects, we adopt the optical clas- 
sifications from IVeron-Cetty & VeronI ( |2OO60 . which include 
three catalogues: 'quasars', 'AGNs' and 'Blazars'. In these 
catalogues, quasars are defined as objects that are starlike or 
with a starlike nucleus, broad emission lines, and brighter 
than absolute magnitude Mb=-23. AGNs include Seyfert Is, 
Seyfert 2s and LINERs that are fainter than Mb=-23. Seyfert 
Is are further divided into five subgroups: Seyfert 1.0, 1.2, 
1.5, 1.8 and 1.9 based on the appearance of the Balmer lines. 
For our study, we compiled a quasar sample including all ob- 
jects from the 'quasar' catalogue that have been classified as 

51, S1.2, SI. 5 or Sin (Narrow Line Seyfert Is) from their 
optical emission line properties; a Seyfert 1 sample from the 
'AGN' catalogue as sources that have optical spectroscopic 
classifications of SI, SI. 2, SI. 5 and Sin; and a Seyfert 2 
sample from the 'AGN' catalogue that have classifications of 

52, Slh (Seyfert 2s with broad Balmer lines in the polarized 
light) and Sli (Seyfert 2s with broad Pa (3 in the near IR). 
Objects that are classified as SI. 8 or SI. 9 have weak broad 
Balmer lines, and their identifications as Seyfert Is depends 
very much on the signal-to-noise of the spectra. Therefore, 
we do not include them in our Seyfert 1 sample. We also do 
not consider objects classified as LINERs in our sample. 

Our ULIRG sample includes all objects observed in pro- 
gram 105 with all four low resolution modules and having 
z < 0.5. It is a loosely defined sample, ch osen primarily on th e 
basis of bolometric infrared luminosity dArmus et al.ll2007h . 
Some ULIRGs also show optical char acteristics of quasars, 
Seyfer t Is or Seyfert 2s according to IVeron-Cetty & VeronI 
yOOo), and these are also included in the three AGN samples. 
The final sample includes 24 quasars (including 3 ULIRGs), 
45 Seyfert Is (including 7 ULIRGs), 47 Seyfert 2s (including 
8 ULIRGs), and 98 ULIRGs, in total 196 separate sources. 

The s pectra were extract ed using the SMART analysis 
package (Higdon et alj |2004l) . Extractions were done differ- 
ently for observations done in staring mode and in mapping 
mode. For staring mode observations, we extract spectra from 
the DROOP products provided by the Spitzer Science Center 
in Pipeline version 11.0.2 and 14.0, and subtract background 
by differencing the SLI and SL2 or LLl and LL2. The spec- 
tra are then calibrated using the IRS standard stars HD17351 1 
(5.2-19.5/im) and ^Dra(19.5-38.5/im). The mapping mode 
spectra are extracted using the default column extraction in 
SMART, utilizing the pipeline BCD products and extract- 
ing only the central image of each observation after subtract- 
ing the outermost image as background (see Buchana n et alj 
(|T006) for technical details). The final spectra were stitched 
and scaled to the LLl spectrum. 

3. RESULTS 

The mid-IR spectra of the objects in the four samples 
show great variety. The average spectra of the four sam- 
ples are presented in Figure [Tf. The average quasar spec- 
trum clearly shows silicate features in emission at both 10 
and 18/im. There are very weak indications of PAH emis- 
sion at 6.2, 7.7 and 11.3/im, and weak emission lines from 
[SIV]10.5/xm, [NeII]12.8Mm, [NeV]14.3/xm, [NeIII]15.5/zm, 
[SIII]18.7/im, [NeV]24.3/im and [OIV]25.9/im. PAH emis- 
sion is also detected in the average spectrum of PG quasars, 
eve n for those that do not s how PAHs in their individual spec- 
tra dSchweitzer et al.|[2006 ). For the Seyfert 1 average spec- 

^ The average spectra are available from http : // ; which will be com- 
pleted after the paper's been accepted 




Fig. 1 . — The average spectra of QSOs (red), Seyfert Is (green), Seyfert 2s 
(blue) and ULIRGs (black). The average spectra are obtained by averaging 
after first nomalizing each spectra at lA.S/im in the rest frame. The 10/xm 
silicate strength measured from the average spectra are: +0.20, -0.21, -0.54, 
and -1.44 for quasars, Seyfert Is, Seyfert 2s and ULIRGs, respectively. 

trum, it is not clear if the 10/im silicate feature is in emis- 
sion or absorption. All the PAH features and atomic emission 
lines are stronger compared with the quasars. The average 
spectrum of Seyfert 2s indicates even stronger PAH emission 
features. The silicate feature at 10/im is clearly in absorption. 
The H2 lines are stronger, but the atomic lines are similar to 
those in the Seyfert Is. Stronger P AHs in Seyfert 2s com - 
pared with Seyfert Is are also seen in lBuchanan et al.l (l2006l) . 
The spectrum of the average ULIRG shows the strongest PAH 
emission compared with the AGNs and also shows clear ab- 
sorption at both the 10 and 18 fim silicate features. The H2 
lines are also more pronounced than those in AGNs. The 
mid-IR spectral slope decreases significantly from ULIRGs 
to quasars. 

To quantify the silicate strength at ^ lO/um, we define the 
silicate strength 5 10 as: 



5io = ln 



/co«f( 10/im)' 



(1) 



where /o;,j(10/im) is the observed flux density at the peak of 
the 10/im feature, and /co,„(10/im) is the continuum flux at the 
peak wavelength, extrapolated above the silicate absorption 
feature (or below the emission feature). For deter mination of 
/£P,„(10/tm), we adopt the methods described in Spoo n et aP 
(|2007). This uses three different techniques for PAH-weak 
spectra, PAH-dominated spectra, and absorption-dominated 
spectra to avoid the contamination of various PAH emission 
features to the continuum and to make best use of uncontam- 
inated continuum regions. 

In Figure 2, we show the distribution of the 10/im sili- 
cate strengths of quasars, Seyfert Is, Seyfert 2s and ULIRGs. 
For objects with silicate strength close to (the bin cross- 
ing Su) = in Figure 2), identification of the feature as emis- 
sion or absorption is strongly affected by the S/N of the spec- 
tra and the uncertainty of the continuum determination. We 
ignore the objects in this bin for further statistical analy- 
sis. The majority of quasars have the 10/im silicate feature 
in emission, but 2 cases clearly show the feature in absorp- 
tion. These include one ULIRG (IRAS 00275-2859) and a 
source from the 2MA SS-selec ted red AGN sample (2MASSi 
J1258074H-232921) (Shi et al. 2006). Their optical spectra 
clearly show broad Balmer emission and blue continuum 
dZhenget all 120021; ISmifli et al.1 l2002h . which do not distin- 
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Fig. 2. — The distribution of tlie lOfim silicate strength, 5io, for QSOs, 
Seyfert Is, Seyfert 2s and ULIRGs. Silicate absorption increases to the right. 
In QSOs, Seyfert Is and Seyfert 2s, the distributions of sources that are also 
in the ULIRG sample are shaded. The averages of the lO/jm silicate strengths 
are 0.20 for quasars, -0.18 for Seyfert Is, -0.61 for Seyfert 2s and -1.56 for 
ULIRGs. The range of the silicate strengths are < 5io < 0.4, -0.7 < Sio < 
0.3, -1.6 < 5io < 0, — 4 < Sio < 0.2 for most quasars, Seyfertl, Seyfert 2s 
and ULIRGs. 

guish them from typical quasars. But the fact that they are 
selected in the IRAS or 2MASS red AGN sample which se- 
lect for dusty sources agrees with their having silicate features 
in absorption. For quasars that have silicate emission, the fea- 
tures are weak. The strongest silicate strength in quasars has 
a value of ^ H-0.75 (PG1351H-640), but the remainder all have 
silicate strength less than 0.4. The average silicate strength in 
the quasar sample is 0.20. 

In the Seyfert 1 sample, even more objects show silicate 
absorption at 10 fim. Excluding objects that have measured 
silicate strength close to 0, the number of objects that have 
silicate in emission is about the same as the number of ob- 
jects that have silicate absorption. All of those in emis- 
sion have strengths less than 0.35, and the average silicate 
strength in the sample is -0.18. There are four Seyfert 1 
sources that show silicate absorption deeper than -1. They 
are 2MASSiJ1307006H-233805, 2MASSiJ1453315H-135358, 
IRAS14348-1447 and UGC5101. However, their opti- 
cal spectra do no t show unambiguou s broad Balmer line s 
(ISmith et alJl200l IVeiUeux et al.lll99l ISanders et alJflQSSh . 
and it is possibly arbitrary to classify them as Seyfert Is. 
The deepest silicate absorption for sources with unambigu- 
ously identified broad Balmer lines is -0.65, from Mrk231 
dWeedman et al. !'20b5). We particularly note that the silicate 
strength in quasars is more dominated by emission than in 
Seyfert 1 AGNs, and that the luminosities of the quasars as 




0.1 

l!5.5Mm)/l(14.5;im) 



Fig. 3. — The coiTelation between the observed strength of the silicate 
feature, Sio, and various IR slope. The open circles are ULIRGs, squares are 
Seyfert Is, crosses are Seyfert 2s and filled circles are quasars. The color 
code is the same as in Figure 1 . 

shown in Figure 4 are systematically greater than for Seyfert 
Is. 

Almost all of the Seyfert 2s in the sample show the 
10 /im silicate feature in absorption except for one case: 
IRAS F01475-0740. The average silicate strength of the sam- 
ple is -0.61. Silicate emis sion is also detect ed in the average 
spectra of type II QSOs ( Sturm et al.l 120061) . although IRAS 
F01475-0740 in our sample is not a highly luminous source. 

The ULIRG sample is dominated by sources with strong sil- 
icate absorption with 10/im silicate strength -4 < Sio < 0.2. 
There are 2 cases of weak silicate emission: Mrkl014 and 
IRAS 07598H-6508, but their optical spectra place the former 
in the quasar and the latter in the Seyfert 1 samples. The av- 
erage of the lOfim strength in the ULIRG sample is -1 .56. 

In Figure 3, we show the comparison between the infrared 
spectral slope and the silicate strength. There is a clear cor- 
relation of the silicate strengths with the IR slopes as mea- 
sured from the flux ratio of 14.5/xm to 27.5/im, a result that 
is also evident in Figure 1. But there is less correlation of 
the silicate strength with the slope measured from the flux 
ratio at 5.5/^m and 14.5/xm or 5.5/im and 27.5/im . This sug- 
gests that the silicate feature arises primarily within the cooler 
dust responsible for the mid-IR emission of deeply absorbed 
ULIRGs, and that the near-IR emission does not arise from 
the same dust. iLevenson et al.l (2007) suggest that absorption 
features as strong as those observed require a nuclear source 
to be deeply embedded in a smooth distribution of dusty mate- 
rial that is both geometrically and optically thick. Given their 
extreme obscuration, the near-IR emission would be heavily 
attenuated. The observed near-IR emission would thus have 
to have a different origin, and not fill in the deep silicate ab- 
sorption trough. 

4. DISCUSSION 

Our sample is not complete in any sense, because it includes 
a variety of objects selected with different criteria. However, 
with 196 objects, this is the largest available dataset of local 
AGNs and ULIRGs observed by the IRS which enables us to 
measure the distribution of the silicate strength. The distri- 
bution of the silicate strengths should be representative of the 
different categories of objects in the local universe. 

The results in Figure 2 are important for understanding 
the faint sources (fy(24/.tm) ~ 1 mJy) at high redshift which 
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Fig. 4. — The correlation between the observed strength of the siHcate 
feature, 5io, and the IR luminosity at 5.5fim . The symbols and color code 
ai'e as in Figure 3; silicate absorption increases to the right. 

have large values of IR/opt for which S/N is poor and fea- 
tures must be strong for a redshift to be assigned. Of the 53 
sources with redshifts derived from silica te absorption as re- 
ported in |H ouck et al. (2005), IWeedman et all (12006b), and 
IWeedmanet al. (2006c) , 44 sources had silicate absorption at 
least as deep as in Mrk 231, for which the silicate strength is 
-0.7 (the absorbed continuum is at a level of 50% of the con- 
tinuum extrapolated above the absorption). 16 of 76 sources 
observed (21%) had no apparent absorption, to a limit for Sm 
of - -0.3. Figure 2 indicates that only 10% of ULIRGS have 
silicate absorption shallower than -0.3 (Sio > -0.3), implying 
that ULIRGS are even more absorbed than the faint sources 
characterised by extreme IR/opt, although the absence of fea- 
tures in the optically faint sources could sometimes arise be- 
cause features are redshifted out of the IRS spectral range. 

Figure 4 shows the 5.5/im luminosity compared with the 
lOfim silicate feature strength for quasars, Seyfert Is, Seyfert 
2s and ULIRGs. It is particularly useful for interpreting the 
sources at high redshift, which must be extremely luminous in 
order to fall within samples having f^(24fim) 1 mJy and z 
^ 2. For example, the most lu minous source yet fo und among 
these samples is source #9 in Houck et al.l (l2005h . which has 
luminosity of 7 x lO'^L© at 5.5/im, and exceeds any source in 
our current sample. All of the sources at z ~ 2 reported so far 
have luminosities above ^ 2 x lO'^L© at rest wavelength of 
~ 6/im. Figure 4 indicates that both quasars and ULIRGS can 
have luminosities that match these high redshift sources. High 
redshift quasars have been discovered in the Spitzer 24/im 



samples, but these quasars have /? < 21 .7 dBrown et al.ll2006h 
and therefore have much brighter optical magnitudes than the 
obscured sources with extreme IR/opt. This means that exam- 
ples at high redshift similar to the "quasars" in Figure 4 have 
been discovered in the infrared, but they are not heavily at- 
tenuated in the optical. This implies that the local analogues 
to the optically faint sources at high redshift would be those 
luminous ULIRGs in Figure 4. 

The nature of the dust distribution in various types of AGN 
and the relation of this to a unification scheme is uncer- 
tain. For example, some silicate emission from AGNs may 
arise on scales la rger than the torus of the unification scheme 
(ISturm et al.ll2005 ). The distributions of the silicate strengths 
presented here are intended to give statistical meaning to un- 
derstanding the dusty structure of AGNs and ULIRGs, rather 
than interpreting individual objects. For example, the aver- 
age spectra provide a practical means to identify fundamen- 
tal differences among the nuclear dusty environments of the 
different families. A significant number of ULIRGs show 
very deep silicate absorption (Sm < -2.5), in strong contrast 
to the quasars of similar luminosity which show no obvi- 
ous sign of absorption and also in contrast to the classical 
optically-classified AGNs. This reflects a different geome- 
try of the dusty mate rial in these deeply absorbed ULIRGs 
dLevenson et al.l l2007). The distribution also reveals a signifi- 
cant number of silicate absorptions in Seyfert Is and a few in 
quasars, and one silicate emission in a Seyfert 2. These results 
are not expected from the simplest unification scenario, al- 
tho ugh several AGN models can accomod ate such exceptions 
(eg.'Nenko va et al.ll2002l:lEfstafliioul2006 . M. Nenkova et al., 
in preparation). The fraction of these exceptional cases and 
the ranges of their silicate strengths provide important con- 
straints on these models. It should be noted that the fraction 
of Seyfert Is having weak silicate absorptions is much larger 
than the fraction of Seyfert 2s with silicate emission, which 
requires a dust geometry that allows the broad-line region to 
be easily visible even when there is significant optical depth 
through the dust clouds themselves. 
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